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Abstract. Let E„ = [xi = 1, x,- + Xj = x^, Xi ■ Xj = Xk ' i, j,k e {\, . . . ,n}}. There 
is an algorithm that for every computable function / : N ^ N returns a positive 
integer m(/), for which a second algorithm accepts on the input / and any integer 
n > m{f), and returns a system S Q E„ such that S has infinitely many integer 
solutions and each integer tuple {xi, . . . , x„) that solves S satisfies xi = f{n). For 
each integer n> 12 we construct a system S Q E„ such that S has infinitely many 

integer solutions and they all belong to Z" \ [-2 ,2 ]". 
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We present a general method for constructing small systems of Diophantine 
equations which have only very large integer solutions. Let 0„ denote the follow- 
ing statement 

Vjci,...,x„ eZ3yi,...,y„ eZ 



A 



(2^' <|xi|^(k|<|yi|V...V|.^i|<|3;„|)) 

(Vz, j,ke{l,...,n}(xi + xj = Xk ^ yi + yj = yk)) A (1) 
Vz, j,ke {l,...,n} {xi ■ xj = Xk ^ yi ■ yj = yt) (2) 
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For n>2, the bound 2^ cannot be decreased because for 

2^ X ,2^ ,...,256,16,4,2) 
the conjunction of statements (1) and (2) guarantees that 

(ji, . . . , j„) = (0, . . . , 0) V (yi, . . . = , 22""', 22"'', . . . , 256, 16, 4, 2) 

The statement VnO„ has powerful consequences for Diophantine equations, 
but is still unproven, see [5J. In particular, it implies that if a Diophantine equation 
has only finitely many solutions in integers (non-negative integers, rationals), then 
their heights are bounded from above by a computable function of the degree and 
the coefficients of the equation. For integer solutions, this conjectural upper bound 
can be computed by applying equation (3) and Lemmas |2] and |7l 

Observation. For all positive integers n, m with n < m, if the statement 0„ 
fails for {x\,. . . ,Xn) 6 Z" and 2 <\xi\ <2 , then the statement <t>m fails for 

{ Xi, . ^. , Xi , X2, . . . , Xn) 6 Z™. 

m-n+l times 

By the Observation, the statement VnO„ is equivalent to the statement 'in^'n, 
where denote the statement 

^xu...,Xn eZ3yu...,yn eZ 

(22" ^ < Ixil = max(|xi|, . . . , |x„|) < 22" =^ {\xy\ < \yi\ V ... V \xi\ < \y„\)) A 

(yi,j,ke {l,...,n} (jc,- + xj = Xk ^ yi+yj = yk)) A 

'^i,j,ke {\,...,n] (xi ■ xj = Xk ^ yt-yj = yk) 

In contradistinction to the statements 0„, each true statement ^„ can be confirmed 
by a brute-force search in a finite amount of time. 

The statement 

Vn Vxi , . . . , x„ e Z3yi, . . . ,y„ e Z 
(2^"'' <\x,\^\x,\<\y,\)A 
(ii,j,ke {\,...,n]{xi + Xj = Xk^ yi+yj = yk)) A 
Vz, J, ^ e {1, . . . , n} {xi ■ Xj = Xk ^ yt ■ yj = yk) 
strengthens the statement VnO„ but is false, as we will show in the Corollary. 
Let 

En = {xi = 1, Xi + Xj = Xk, Xi • Xj = Xk '• i, j,k G { 1, . . . , n}} 
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To each system S Q E„ we assign the system S defined by 

{S\{Xi = 1 : /6 {l,...,n}})U 
{Xi ■ Xj = Xj : i, j e {I, . . . ,n} and the equation x,- = 1 belongs to S } 

In other words, in order to obtain S we remove from S each equation Xi = 1 and 
replace it by the following n equations: 

Xj • Xi = Xi 

Lemma 1. For each system S Q E„ 

{(xi , . . . , 6 Z" : {xi, . . . , x„) solves S } = 
{(xi, Xn) eZ" : (jci, . . . , x„) solves S } U {(0, . . . , 0)} 

Lemma 2. The statement 0„ can be equivalently stated thus: if a system S Q E„ 
has only finitely many solutions in integers xi,...,x„, then each such solution 

(xi, . . ., Xn) satisfies \xi\, . . ., \x„\ < 2 

Proof. It follows from Lemma [B □ 

Nevertheless, for each integer n> 12 there exists a system S Q E„ which has 

infinitely many integer solutions and they all belong to Z" \ [-2 ,2 ]". We 
will prove it in Theorem [1] First we need a few lemmas. 

Lemma 3. If a positive integer n is odd and a pair {x, y) of positive integers solves 
the negative Pell equation x^ - dy^ = -1, then the pair 




[ 2 • 2V3 J 

consists of positive integers and solves the equation x^ - dy^ = -1. 

Lemma 4. pp. 201-202, Theorem 106]) In the domain of positive integers, 
all solutions to x^ - 5y^ = - 1 are given by 

(2+ Vs)^^^^ =;c + >'V5 
where k is a non-negative integer. 



3 



Lemma 5. The pair (2, 1) solves the equation - 5y^ = - 1. If a pair {x, y) solves 
the equation x^ - 5y^ = -1, then the pair {9x + 20j, Ax + 9y) solves this equation 
too. 

Lemma 6. Lemma \5\ allows us to compute all positive integer solutions to 
x^-5y^ = -1. 

Proof. It follows from Lemma 51 Indeed, if ^2 + Vs) = x + y Vs, then 
(2 + ^^5f^^ = (2 + ^^5f ■ (l + Vs)^*"^^ = 
(9 + 4^^5)■(x + y^^5) = (9x + 20y) + {4x + 9y) Vs 

□ 

Theorem 1. For each integer n > 12 there exists a system S Q E„ such that S 
has infinitely many integer solutions and they all belong to Z" \ [-2 , 2 ]". 

Proof. By Lemmas |4]-[6l the equation u^ - 5v^ = -1 has infinitely many solutions 
in positive integers and all these solutions can be simply computed. For a positive 
integer n, let iu{n), v{n)) denote the n-th solution to u^ - 5v^ = -1. We define S as 

^1 = 1 X\ -\- X} = X2 X2 + X2 — X^ X\ + JC3 = JC4 

Xi\ • X4 = ■ X5 = X(, X(y ■ X-] = JCg • JCg = X() 

X\0 ' — -^11 + .^1 = X4 • Xc) = X\2 

Xn ■ ^12 = ■ ^13 = ■ ■ ■ ^n-\ " ^n-\ = 

The first 1 1 equations of S equivalently expresses that x\q - 5 ■ x\ = -I and 625 
divides xg. The equation x\q - 5'^ ■ x^ = -I expresses the same fact. Execution of 
the following MuPAD code 

x:=2: 
y:=l: 

for n from 2 to 313 do 

u:=9*x+20'-y: 

v:=4*x+9*y: 

if igcd(v, 625)=625 then print(n) end_if: 

x:=u: 

y:=v: 
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end_£or : 
£loat(u^2+l) ; 
float(2^(2^(12-l))) ; 

returns only n = 313. Therefore, in the domain of positive integers, the minimal 
solution to x]q - 5"^ ■ Xj = -1 is given by the pair ^jcio = m(313), xj = ^^^25^ j- 
Hence, if an integer tuple (xi, . . ., Xn) solves S , then Ixg] > v(313) and 

xn = + 1 > uOnf + 1 > 2^^^ ^ 

The final inequality comes from the execution of the last two instructions 
of the code, as they display the numbers 1.263545677e783 and 

3.231700607e616. Applying induction, we get x„ > 2 . By Lemma |3J 
(or by [|6l p. 58, Theorem 1.3.6]), the equation x\q - S'^ ■ x^ = -\ has infinitely 
many integer solutions. This conclusion transfers to the system S . 

□ 

J. C. Lagarias studied the equation x^ - dy^ = -\ for d = 5^"^^ , where 

n = 0, 1, 2, 3, His theorem says that for these values of d, the least integer 

solution grows exponentially with d, see [2, Appendix A]. 

The next theorem generalizes Theorem [1] But first we need Lemma|7]together 
with introductory matter. 

Let D{x\ , . . . , Xp) e Z[xi , . . . ,Xp]. For the Diophantine equation 
2 ■ D{xi, . . . ,Xp) = 0, let M denote the maximum of the absolute val- 
ues of its coefficients. Let T denote the family of all polynomials 
W(xi, . . . ,Xp) 6 Z[xi, . . . ,Xp] whose all coeflicients belong to the interval 
[-M,M] and deg(W, x,) < dt = deg(D, x,) for each / 6 {l,...,p}. Here we 
consider the degrees of W(xi , Xp) and D{xi , Xp) with respect to the 
variable Xj. It is easy to check that 

card(r) = (2M + l)*^^i + 1) • • • • • (^^p + 1) (3) 

We choose any bijection t : {p + I, . . . , card(7~)} — > T \ {xi, . . . , Xp}. Let *K 
denote the family of all equations of the form 

Xi = \, Xi + Xj = xt, Xi ■ Xj = Xk {i,j,ke{l,..., card(7~)}) 

which are polynomial identities in Z[xi, . . . , Xp] if 

V5 6 {p + 1, . . . , card(7~)} x, = t(s) 
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There is a unique q £ {p + I, . . ., card(7~)} such that T(q) = 2 • D(xi, Xp). For 
each ring K extending Z the system implies 2 • D(xi , Xp) = Xq. To see this, 
we observe that there exist pairwise distinct to,...,tm&T such that m> p and 

?0 = 1 A ?i = Xi A ... A tp = Xp A t„ = 2 • D(xi ,.. ., Xp) A 

'ii e{p + l,...,m]3j,ke{0,...,i-l} itj + tk = tiV u + h = tj V tj • tk = td 

For each ring K extending Z and for each xi, . . . , Xp £ K there exists a unique tu- 
ple (Xp+u XcaidiT)) e ^':^d(r)-p gyj,]^ jj^^j jj^g j^pig _ ^ Xp, Xp+u Xc^diT)) 

solves the system "H. The sought elements Xp+i, Xcardcr) are given by the for- 
mula 

V5 G {p -I- 1, . . . , card(7~)} = t(s)(xi, ...,Xp) 

Lemma 7. 77je system H U {Xg + Xg = x^} con &e simply computed. For each 
ring K extending Z, the equation D(xi, . . .,Xp) = is equivalent to the system 
9i U {xq + Xq = Xq) Q £'card(T')- Formally, this equivalence can be written as 

Vxi, ...,Xp£K {p{Xu ...,Xp)=Q <^ 3Xp+i, XcardiT) € K 

(xi, . . . , Xp, Xp+i, Xcard(7-)) solves the system 'Hu{xq + Xq = x^}) 

For each ring K extending Z and for each xi,...,Xp £ K with D(xi, . . . , Xp) = 
there exists a unique tuple fx^+i, . . . ,Xcard(r)) ^ ^^''^'^^"^ such that the tuple 
(xi, . . . , Xp, Xp+i, . . . , Xcard(r)) solvcs the system 9^ U {Xq + Xq = Xg}. Hence, for 
each ring K extending Z the equation D{xi , . . ., Xp) = has the same number of 
solutions as the system 9i U {xg + Xg = Xg}. 

Putting M = M/2 we obtain new families T and H. There is a unique q £ 
{l,...,card(7")} such that 

{q£ {l,...,p] A Xg = D(xi,...,Xp)) V 

{q£{p + l,..., card(7~)} A T(q) = D(xu Xp)) 

The new system U {x^ -i- x^ = x^} is equivalent to D(xi ,...,Xp) = and can be 
simply computed. 
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The Davis-Putnam-Robinson-Matiyasevich theorem states that every recur- 
sively enumerable set A1 c N" has a Diophantine representation, that is 



(ai, . . . , a„) 6 A1 <=^ 3^1, . . .,Xm eN W(ai, . . . ,an, xi, . . . , x,n) = 

for some polynomial W with integer coefficients, see f3\ and (\]. The polyno- 
mial W can be computed, if we know a Turing machine M such that, for all 
(ai, . . . , a„) 6 N", M halts on (ai, . . . , a„) if and only if (ai, . . . , a„) G M, see Q 
andffl. 

Theorem 2. There is an algorithm that for every computable function / : N ^ N 
returns a positive integer m{f), for which a second algorithm accepts on the 
input f and any integer n > m(f), and returns a system S Q En such that S has 
infinitely many integer solutions and each integer tuple (xi, . . . , Xn) that solves S 
satisfies Xi = f{n). 

Proof. By the Davis-Putnam-Robinson-Matiyasevich theorem, the function / 
has a Diophantine representation. It means that there is a polynomial 
W{x\ , JC2, . . . , Xr) with integer coefficients such that for each non-negative inte- 
gers XuX2, 



By the equivalence (El) and Lagrange's four-square theorem, for each integers 
X\, X2, the conjunction {x2 > 0) A {xi = f{x2)) holds true if and only if there exist 
integers a, b, c, d, a,j3, y, 6, x^, x^j, ^3 2, x^^^, X3 4, . . . ,Xr, x^j, Xr^, Xr^3, Xr^ such that 



W^ixi, X2, X3,..., Xr) + (xi -a^-b^-c^- d^f + {x2 -a^-/3^-y^- 6^f+ 



By Lemma|71 there is an integer s > 3 such that for each integers xi, X2, 



where the formula ^{xi,X2, x^,, . . ., x,) is algorithmically determined as a conjunc- 
tion of formulae of the form Xi = 1, Xi + Xj = xt, Xi ■ Xj = Xk (i, j,k e {\, . . . , s}). 
Let m(/) = S + 2s, and let [•] denote the integer part function. For each integer 
n > m{f). 



= f(X2) 



3^3, . . .,Xr & M W(Xi,X2, X3, . . ., X,.) = 



(El) 





-J - 4 - ^ > m(/) - 



m(/) 



-4- s> m(f) - 



m(f) 



-A-s=0 



n - 



2 



2 
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Let S denote the following system 



all equations occurring in JCa, ...,Xs) 
n - 1^1 j - 4 - s equations of the form z, = 1 

h 

ti + h 



1 



< 



^[2]-i + h 



w 

y (if n is even) 
1 (if n is odd) 



w + y 
y + y 



y 

u + u 



V 



with n variables. By the equivalence (E2) , the system S is consistent over Z. The 
equation u + u = v guarantees that S has infinitely many integer solutions. If an in- 
teger n-tuple (xi,X2, ... ,Xs, ... ,w,y,u,v) solves S , then by the equivalence (E2 ) , 



Corollary. There is an algorithm that for every computable function / : N ^ N 
returns a positive integer m(f), for which a second algorithm accepts on the 
input f and any integer n > m(f), and returns an integer tuple (xi, . . . , x„) for 
which xi = f(n) and 

(4) for each integers yi,...,yn the conjunction 



'ii,j,ke{l,...,n} (xi ■ xj = Xk ^ yt • yj = yu) 
implies that x\ = yy. 

Proof. Let <„ denote the order on Z" which ranks the tuples ( jci , . . . , x„) first 
according to max(|xi|, . . . , |x„|) and then lexicographically. The ordered set 



Xi = 



f(x2) = /(w + y) = / (2 • [^] + = f(n) 



□ 



(Vie{l,...,n} (x, = 1 ^3;, = 1)) A 
(v/, j,kG {l,...,n} (xi + Xj = Xk^ yi + yj = yk)) 



A 
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(Z", <„) is isomorphic to (N, <). To find an integer tuple (xi, . . . , jc„), we solve 
the system S by performing the brute-force search in the order <„. 

□ 

If n > 2, then the tuple 

2^ ,2^ ,.. .,256,16,4,2,1] 

has property (4). Unfortunately, we do not know any explicitly given integers 
xi, . . .,Xn with property (4) and \xi\ > 1 
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